Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 24 March 2010 (MN M£X style file v2.2) 



Fundamental parameters of RR Lyrae stars from 
multicolour photometry and Kurucz atmospheric models. 
I. Theory and practical implementation 
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ABSTRACT 

A photometric calibration of Kurucz static model atmospheres is used to obtain pa- 
rameters of RR Lyrae stars: variation of stellar angular radius 1?, effective temperature 
T G , and gravity g e as a function of phase, interstellar reddening E(B — V) toward the 
star, and atmospheric metallicity M. Photometric and hydrodynamic conditions are 
given to find the phases of pulsation when the quasi-static atmosphere approximation 
(QSAA) can be applied. The QSAA is generalized to a non-uniformly moving spher- 
ical atmosphere, and the distance d, mass M, and atmospheric motion are derived 
from the laws of mass and momentum conservation. To demonstrate the efficiency of 
the method, the UBV(RI)c photometry of SU Dra was used to derive the following 
parameters: [M] = -1.60 ± .10 dex, E(B - V) = 0.015 ± .010, d = 663 ± 67 pc, 
M = (0.68 ± .03)A4 Q , equilibrium luminosity, L cq = 45.9i9.3Lo, T cq = 6813±20 K. 

Key words: stars: variables: RR Lyr - stars: fundamental parameters - stars: oscil- 
lations - hydrodynamics - stars: individual: SU Dra - stars: atmospheres 



1 INTRODUCTION 

Although the determination of the fundamental parameters 
of RR Lyrae (RR) stars is interesting in itself, it is also im- 
portant from a practical point of view, because RR stars play 
a considerable role in establishing galactic and extragalactic 
distance scales. The Preston index and spectroscopic obser- 
vations are used to determine their atmospheric metallicity 
[M] and interstellar reddening E(B — V). Since confirmed 
RR type components are not known in binary systems, the 
mass determination is based on both stellar evolution and 
pulsation theories. Due to the uncertainty of parallax data, 
mostly the Baade-Wess elink (BW) method is used to infer 
their distance d (|Smithll"l995l ). In the BW analysis, and to 
determine [M] and E(B — V), the quasi-static atmosphere 
approximation (QSAA) is employed to interpret photometry 

and spectroscopy. 

The QSAA was introduced by iLedoux fc Whitney! 

"The simplest approach is to assume that at each 
phase, the atmosphere adjusts itself practically instanta- 
neously to the radiative flux coming from the interior and 
to the effective gravity g c 



g c = GM.BT + R 



(D 



where R, and R are the instantaneous values of the radius 
and acceleration, which is supposed uniform throughout the 
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atmosphere", G is the Newtonian gravitational constant, M. 
is the stellar mass, and dot is a differentiation with respect 
to time t. "One may then build a series of static model atmo- 
spheres, " and select one of them at each phase by spectro- 
scopic or photometric observations. Its flux, colours, effec- 
tive temperature T c , and surface gravity g c are accepted as 
the atmospheric parameters of that phase providing basis 
for determination of other parameters like angular radius, 
mass, distance, etc. 

The subject of the present paper is the QSAA and its 
generalization to a non-uniform atmosphere. We will inves- 
tigate the QSAA from the point of view of atmospheric 
emergent flux and hydrodynamics. By comparing the ob- 
served colour indices with those of s tatic model at mospheres 
l|Castelli. Gratton fc Kurucj|l997l. lKuruczlll997l ) we select 
the phases when they coincide. Considering hydrodynamics, 
we do not construct a consistent dynamic model of an RR 
atmosphere. However, we find a better description of the 
pulsating atmosphere if we characterize it by pressure and 
density stratifications in addition to the two parameters R 
and R. We determine the fundamental parameters M and 
d from the hydrodynamic considerations without using the 
BW method or theories of stellar evolution and pulsation 
at all. Our method uses photometry as observational in- 
put; spectroscopy and radial velocity observations are not 
needed. 

In Section 2, conditions of QSAA are formulated for a 
spherically pulsating compressible stellar atmosphere with 
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velocity gradient. Practical methods are described to deter- 
mine [M], E(B - V), T e (<p), and log fle (¥>) from UBV(RI) C 
colours of Kurucz atmospheric models, ip being the pulsa- 
tion phase. The laws of mass and momentum conservation 
are used to determine the mass and distance of the pulsat- 
ing star. Section 3 presents the results obtained from the 
UBV(RI)c photometry of SU Dra. Discussion and conclu- 
sions are given in Sections 4 and 5. 



2 THE QUASI-STATIC ATMOSPHERE 
APPROXIMATION 

First, we describe a photometric method to select the phases 
in which QSAA can be regarded as a good approximation 
from a point of view of fluxes. Next, we introduce the time- 
dependent pressure and density stratifications of the selected 
static models in the laws of mass and momentum conserva- 
tion to find the phases in which QSAA can be regarded as 
a more or less good approximation from a hydrodynamic 
point of view as well. Using the phases of valid QSAA from 
photometric and hydrodynamic points of view, we determine 
mass, distance, and infer the internal motions of the atmo- 
sphere with respect to the stellar radius R. In the present 
paper, we define R as where the optical depth is zero, since 
the emergent fluxes of a theoretical model are given for r = 
in any photometric band. Of course the zero boundary con- 
dition is a mathematical idealization (|Ledoux fe Whitney! 
Il96lf ). In the following, we tacitly assume that the density 
is not zero at the boundary, i.e. g(R) >0if0<r-Cl. 

2.1 QSAA from photometric point of view 

The functions T e [CIi, Ch, [M], E(B - V)] and 
ff e [CIi,CI 2 , [M],E(B - V)] of the Kurucz atmospheric 
models form a suitable grid for interpolation in the ranges 
6000 <T C < 8000 and 1.5 < log g c < 4.5 if the colour indices 
CIi,Cl2 are selected appropriately from UBV(RI)c pho- 
tometry. For fixed [M] and E(B — V) the intersection of the 
two functions {T c (l) (log 5c , CIi, CI 2 , [M], E{B - V))} i=lj2 
as a function of logg c gives a pair of T D (ip), log g c (f) for a 
given pair of Cli(p), Cl2(p)- Although it is more practical 
to us e U - IB + V instead of U - B l|Barcza fc Benkol 
2009), the other common colour indices, e.g. U — V etc., 
will also be kept. 

Using UBV(RI)c photometry, we can construct four 
independent colour indices CI;=i,2,3,4, i.e. (^j = 6 combina- 
tions of colour indices resulting in six pairs of T c (<p), log g e (<p) 
for a given ip. Since the U band covers the Balmer jump, the 
indicator of g c (f), the combinations must contain at least 
one colour index containing U photometry. (^J — — 30 
such combinations can be constructed. The average of the 
30 pairs of T B (tp), log g c {<p) will be accepted for a phase <p. 

Condition I. // the scatters AT a (ip), A log g c (p) of 
T e (<p), log g e (<p) are compatible with the expected scatter from 
the error of the colour indices, QSAA provides a good ap- 
proximation in this phase from the photometric point of 
view. 

The bolometric corrections and UBV(RI)c physical 
fluxes of a model with gi ven T e , logfle, [M], E(B — V) are 
available in tabular form (|Kuruczlll997t ) and they allow to 
determine the angular radius of the star defined by 



0foO = R(v)/d. 



(2) 



Technical details of finding the appropr iate model an d de- 
termining i) at a phas e tp are given in iBarczal l|2003h and 
iBarcza fc Benkol (|2009h Fl 

On the scale provided by the Kurucz atmospheric 
models, the above procedure offers a possibility to de- 
termine reddening and atmospheric metallicity from pho- 
tometry without using spectroscopy. In the shock free 
phases the QSAA is expected to reproduce the colours 
well, therefore, AT^Alogpe must be minimal as a func- 
tion of the (^-independe nt model parameters E(B — V), [M] 
l|Barcza fc B cnko 20C||). Of course the method can also be 
applied for non- variable stars to determine E(B — V), [M], 
T c , log g, tf. 



2.2 QSAA from hydrodynamic point of view 

In the original form of the QSAA (|Ledoux fc Whitnevll96lh , 
the atmospheric motion is described by the parameters 
R(ip) , R(p) , R(p) and the outward or inward accelerations 
are driven by the variable g e (<p) — g s [R(ip)] in (fJJ where 
g B = GM/r 2 . This is essentially a uniform atmosphere ap- 
proximation (UAA) . However, the atmosphere of an RR star 
is a dilute, compressible gaseous system with variable tem- 
perature. Therefore, the introduction of additional param- 
eters containing T c {tp) promises a better description in the 
frame of the QSAA. 

Conside ring the geometry of radial pulsation modes 
(|Smithlll995h and neglecting rotation will result in all com- 
ponents of the vector field velocity v(r),t) being zero except 
for the radial component v(r, t). The viscosity is negligible at 
the velocities occurring in an RR atmosphere. Therefore, the 
momentum conservation is expressed by the Euler equation 
of hydrodynamics: 

dv , dv r \ , 1 d P (tang) / ,n n /on 

where t = Pip, P is the pulsation period. 

The contribution of the neglected tangential motions is 
the average a (tang) = (47T)" 1 dd sin 9 dcf>[v e (dvr/d e) + 
W0 sin" 1 9{dvr/d4>)-vl-vl]/r (|Landau fc Lifehitzlll98Ch . In 
perfect spherical symmetry a^ 118 ^ = 0. The actual pressure 
and density stratifications were divided in two parts p(r, t) = 
p(r,t) +p ( - dyn \r,t) and g(r,t) = g{r,t) + g <dyn) (r,t), where 
p(r, i), g(r, t) denote the values of a static model atmosphere, 
while p( dyn \r, t), g^ dyn '(r,t) are the dynamical corrections 
due to pulsation. 

From the hydrodynamic point of view, the essence of 
the QSAA is that at a given t 



g or 



(4) 



i.e. the r-independent effective gravity of a static model 
atmosphere is introduced in © and p( dyn ) (V, t) = 0, 
p( d yn) _ o. In what follows we drop this latter restric- 
tion and define a dynamical correction term of acceleration 

(dyn) _ 1 dp I dp 



g dr g dr 



(5) 



1 A program package to perform the conversion of colour indices 
and magnitudes to T c , logg. $ is available upon request. 
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to account for the difference of a static and dynamical atmo- 
sphere. Because of the lack of dynamical model atmospheres, 
a ( d y n ) mus t be determined empirically from the quantities 
derived from photometry. 

Eq. © is transformed to the definition equation {TJ 
of QSAA if dv/dr = and a (tang) = 0, a (dyn) = 0, be- 
cause these simplifications imply dv/dt = R and g s (R) = 
GMR~ 2 , i.e. the UAA is valid. As a dynamical equa- 
tion, this simplified equation of motio n is the basis of 
mass determination in BW analyses (e.g. lLiu fe Jane sl ll990l . 
ICacciari. Clementini fc Buscr 1989, etc). 

If the sound velocity is assumed as an upper limit for vg 
and V(p , a*-* 3,115 ' is some 0.1ms -2 . It is negligible in compar- 
ison with the other components of the acceleration at any 
r < R. Typical values are g B < 10ms~ 2 during the pulsation 
cycle of an RR star while g c can exceed 100ms~ 2 when the 
atmosphere is in the state of maximal compression. Thus, 
a (tan g ) w j n k e ne gi ec t e d and 



g(r,t) = g c (t) - g s (r) 



(6) 



is a periodic function of t with zero points. 

Neglecting the r-dependence of the temperature T, 
assuming a constant density upper boundary condition 
g[R(t)], an integration of @ over the interval [r,R] with 
the equation of state of a perfect gas gives the approximate 
density stratification of the model atmosphere: 



g(r, t) = q(R) exp{-/i (R, t) [r - R(t)]} 



(7) 



where ho(R,t) = fig e (t)/lZT(R), 1Z, /i are the reciprocal 
barometric scale height at R, the universal gas constant, 
and average molecular mass, respectively. Essentially we use 
the static model atmospheres to measure — g~ 1 (dp/dr), the 
atmospheric pulsation is driven by g(r,i) — a fdyn \ and the 
thermal processes are represented by the variable ho- 

Condition II. If dv/dt + vdv/dr g, 
i.e. \a fdyn \r, t)\ < 2Ag, QSAA provides a good approx- 
imation in this phase from the hydrodynamic point of 
view. 

This condition formulates that the dynamical excess of 
acceleration in the upper photosphere is smaller than the Ag 
error of g and the same error is assumed for dv/dt + vdv/dr. 
The combination of (JSJ and allows to estimate the 
constant and 0(d) terms of tr dyn . However, the satisfac- 
tion of Condition II can be checked afterwards when the 
terms with other powers of d were determined. The quo- 
tient q(r,t) — a' dyn '/g characterizes the degree of excellence 
of the QSAA. The QSAA is exact from the hydrodynamic 
point of view if q — 0. 

The perfect spherical sym metry means that in Eulcr 
picture (|Pringle fc Kind 120071 1 we have to use ©-([71) with 
a (tan g ) ^ ^ _ q an£ j magg conservation law 



do do dv 2 

■£ + v-^ + £>tt- + -vg = 

dt dr dr r 



(8) 



must be taken into account. With assumption (0 the ana- 
lytic solution of JS} is 

v(r,t) = —air + a + a-ir^ 1 + a-2r~ 2 (9) 

where a\ = Ch^ 1 (dho/dt), ao = R + Ca\(R — 3hg x ), a_i = 
2Ca ho 1 , a_ 2 = Ca-i/i^ 1 , C = 1. C was introduced to 
incorporate UAA by taking a (dyn) (R, t) = and C = 0. The 
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Figure 1. T c , logg c , ho(R) of SU Dra as a function of phase. 
Filled circles indicate the phases when Condition I is satisfied. 



main term in v is R and the convergence of ([9]) is excellent 
because rho S> 1. 

After differentiations of v, ho, and taking r = R, the 
following form is convenient for a numerical solution: 



M(d,ti) - M{d,t 2 ) = 
where M(d, t) 

,(dyn) 



(10) 



[g c (t) - (dv/dt) - v(dv/dr) 



a {dyn> (R,t)]R z /G. It can be solved if there exist two or 
more t intervals satisfying Condition I and a < - dyn ' 1 (R, t) is 
identical for them. The assumed differences a^ dyn \R, ti) — 
a ( dyn \R,t2) ~ and the satisfaction of Condition II must 
be verified afterwards. Eq. (|10p must be solved for different 
phase pairs ((pi, ^2) and the values d and M must be omit- 
ted from the final averaging if a <dyn) (R, <px) « a <dyn) (R, (p 2 ) 
is violated. 

To account for the difference of the local and effective 
temperature, the boundary temperat ure T(R) = 2~ °' 25 T e of 
a gray atmospheric model at r = (Miha las l ll978h will be 
used in ho(R,t). Hydrogen and helium are mainly neutral 
at the boundary temperature of RR stars. Consequently, 
assuming a typical chemical composition means that fj, = 1.3 
is the appropriate choice in ho(R,i). 



3 RESULTS FOR SU DRA 



iBarczal (|2002l 'l derived good quality light curves of the RRab 
star SU Dra by homogenizing 45 years of U BV(RI)c ob- 
servations (see his Table 7). The method described in the 
previous section will now be applied. Instead of smoothing 
the light curves, binned data will be used to keep the results 
as close as possible to observed quantities. The 30 colour 
index pairs containing at least one U were obtained from 
CIi|gi = {U - 2B + V, U - V, U - R C , U - I c , B - V, B - 
Rc,B-I c ,V- Rc, V - Ic, Rc -/cl- Taking [M] = -1.60, 
E(B - V) = 0.015 (|Liu fc Janeslll990l ) the basic quantities 
for (110[) are summarized in Table [T] for different phases. Fig. 
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Table 1. The basic quantities for IllOU at various phases. The 
units are cms _2 ,K,rad X 10 10 ,cm -1 X 10 10 . 
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Figure 2. Upper panels: histogram of 30 logg e , T e values in 
bins of size 0.1, 10K, respectively. Dotted: SU Dra at tp = 0.5, solid 
line: SU Dra at tp = 0.98, dashed: BD +67 708. Lower panels: 
standard errors A log g c and AT C from 30 possible combinations 
of Cli as a function of phase. At phases below the dotted lines 
the scatter can be attributed to observational errors of the colour 
indices. 



1 is a plot of T c (tp), #(<£>), lagg e (<p), ho(R,tp) for one whole 
pulsation. 

AT a (tp), A log g e (tp) are plotted in the lower panels of 
Fig. [2] Assuming random errors of ±0.02 for the colour in- 
dices will result in AT C = ±10K and Alogpc = ±0.04. 
These values are indicated by the dotted horizontal lines. 
Condition I is satisfied in the phase points lying below or 
close to the dotted lines. At phases lying above the dot- 
ted lines, the monochromatic flux of static models and SU 
Dra differs significantly on a level which has noticeable ef- 
fect on the broad band colours U BV(RI)c ■ The assumed 
[M] = —1. 6, E(B — V) = 0.015 w ere verified by a variation 
procedure |Barcza fc Ben ko 2009). In the shock free phases 
p = 0.15,0.5,0.55, minimization of AT e ((p), A log g e (tp) re- 
sulted in [M] = -1.60 ± 0.10, E(B - V) = 0.015 ± 0.01. 

To demonstrate the difference between good and poor 
QSAA, the histograms of the 30 loggr e ,Te values are plotted 
in the upper panels of Fig.[5]for tp = 0.5, 0.98 of SU Dra and 
BD ±67 708. They show normal distributions with small 
scatter for the non-variable BD ±67 708 and SU Dra at 
tp = 0.5. At tp — 0.98 the distribution is almost uniform. 
At the next phase point tp — 1 merely 14 intersections of 
{T e (i) (logffe,CIi,CI 2 ,[M],£(B - V))}i=i,2, i.e. 14 pairs of 
log<; c , T c were found instead of 30 pairs. Thus, at tp « 1 the 
observed colours d iffer s ignificantly from those of any static 
model of iKurucj 1 19971 ). However, it is interesting to note 
that the small errors Alog<; c = 0.03, AT e — 12K indicate 
a phase island at tp = 0.93 with observed and static model 
colours in agreement for some 0.01P w 10 minutes. This 
happens to be the phase of the hump on the light curve when 
the in ward and ou tward motions encounter and produce a 
shock jSmithll 19951) . 
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Satisfaction of Conditions I and/or II is indicated in the last col- 
umn by I and/or II, respectively. 

In the interval 0.92 < tp < 1.05 the atmosphere is in 
a state of maximal compression by the shock coming from 
the sub-photospheric layers and R is nearly minimal. This 
is the rising branch and the start of the descending branch 
in the light curve. Therefore, the values of logg ,T ,i? ob- 
tained from QSSA, if they can be found at all, must be 
considered as a first approximation only. This is reflected in 
large A\ogg e , AT e except for tp w 0.93. 

To get v(r,t) etc for (JTTJJ , #(t) and ho{R,t) were differ- 
entiated by midpoint formulae. Fig. [3^l is a plot of the func- 
tions M(d, t) for the phases t = <pP, tp = 0.15-0.35, 0.5, 0.55 
with a (dyn) (7?,i) = 0. The average and standard error of 
M,d are given in Table [2] from the pairs (tpi = 0.25 and 
p 2 = 0.15,0.3,0.35,0.55), (tpi = 0.35 and <p 2 = 0.55) in 
(|10|) as our best values denoted by [*]. At tp — 0.35 Con- 
dition I is moderately violated but Condition II is satisfied 
and a (dyn) (R,t)/g s (R,t) sa -0.07, therefore, this phase was 
included in getting d, M of [*] . Condition I is satisfied at 
p = 0.2, 0.5, however, these phases had to be excluded from 
the mass and distance determination because of the large 

„(dyn) 



(R, tp = 0.2, 0.5) = 5.8, -3.3 ms^ and q 



-0.52,38, 



respectively. 

Using our best solution for M and d, the radius vari- 
ation, velocities, and the components of acceleration were 
computed in physical units and are plotted in Fig. [3]j-d. Ve- 
locities and accelerations are plotted only for the phases of 
more or less good QSAA (tp = 0.15-0.3,0.5,0.55), including 
the slightly shocked phases tp = 0.35-0.45. 

At the phases tp = 0.15,0.25,0.3,0.35,0.55 



,(dyn) 



-79,8,-16,65,93 cms" 



and Condition 

(dyn) / 



I is satisfied, \a {dyn) (R,t)/g s {R,t)\ < 0.13, i.e. a {dyn > (R, tp) 
is small in comparison with the other acceleration terms 
in ©. |g| w 0.1 is expected from Alog(; c = 0.04, 
q(R,tp) = 0.39, 0.02, 0.06, -0.10, -0.13 were found. 
Thus, Condition II is indeed satisfied. The outlier value 
0.39 is produced by cancellation of dv/dt w —30.5 ms -2 
and vdv/dr ~ ±28.2 ms -2 , i.e. the free fall of the upper 
atmosphere is for a short time almost uniform with constant 
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Figure 3. Panel (a): The functions M(d,cp) in units of solar 
mass and pc, respectively. Dotted, from left to right: cp = 0.2, 0.5, 
dashed: cp = 0.25, lines: from left to right cp = 0.55, 0.35, 0.3, 0.15. 
Panel (b): crosses: R(cp), squares: R(cp) — l/ho(R), in so- 
lar units. Panel (c): crosses: R(cp), squares: v(cp), triangles: 
kg 1 (dv I 'dr)\ r —n. Panel (d): crosses: R(cp), squares: dv/dt + 
vdv/dr, circles: a( d y n ) (R, ip), thick line: g s (R). In Panels (b)-(d) 
the angular variables were converted to absolute radius etc. by 
d = 663 pc and M = 0.68Mq. 



Table 2. Distance and mass of SU Dra from (1 1 1) . 
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C = 0, UAA 



[*] indicates the values which we accept as the best ones. 

velocity and vanishing acceleration —g(R,cp) ~ —2.3 ms -2 . 
Therefore, the density stratification differs from (O and 
this is manifested in the large q. However, this phase could 
also be used to get the best d and M in Table [2] because 
|a (dyn) | < 1ms" 2 . 

In the interval 0.92 < cp < 0.95 the acceleration 
is dominated by g e ~ 89, 112, 148, 102 ms -2 2> g B and 
the atmosphere is practically in standstill. The values 
q = 1.06,-1.68,-0.12,1.44 and a (dyn) (7?, t)/g s (R, t) « 
5,28,-2,6 indicate a (dyn) « 41, 254, -19, 55 ms' 2 and 
| a (dyn)| ^ ffs if g s i s neglected by setting M{d,t) = and 
the data of Table[T]are used with a' dyn ' = 0, the upper limit 
d < 642pc is obtained from cp — 0.94. If this estimation is im- 
proved by a (iyn) (R,cp = 0.94) = -19ms" 2 obtained for the 
distance and mass [*], the upper limit shifts to d < 724pc. 
From the functions .M(724pc,<£ = 0.15,0.25-0.3,0.55) the 
upper mass limit is Ai < 0.84 ± .03 A4q. These upper lim- 
its are in accordance with the values [*]. However, they 
should be treated as indicative only, because they are sensi- 
tive to the actual value of #, i9, dho/dt, d 2 ho/dt 2 . These val- 



ues could be determined by numerical differentiations from 
ho in 0.91 < cp < 0.96 when the QSAA is not optimal 
to obtain all of them. Hence the phase island cp — 0.93-0.94 
could not be included in the mass and distance determina- 
tion because a very transient good QSAA is embedded in a 
strongly shocked interval and the mass M appearing in g s 
plays a secondary role in the dynamics. 

For the sake of completeness it must be mentioned that 
the pair (epi — 0.15 and cp 2 = 0.3) in (|10|l gives d = 1485 pc. 
Adding this value to those involved in [*], d = (800T148) pc 
and M — (1.15 ± .05)^4© is obtained. This distance is 
closer to d — 900 pc su ggested by the H ipparcos paral- 
lax 7r = 1.11 ± 1.15 mas l|Perrvmanlll997f ). However, this 
1485 pc was omitted because it is an outlier point above 3 
sigma level, it exceeds the upper limit 724 pc, and the inac- 
curate knowledge of a (dyn) (fl, 0.15) -a (dyn) (J?,0.3) makes it 
uncertain. 

Close to the shocked phases (i.e. for cp — 0.1, 0.4, 

0. 45, 0.7-0.76, 0.84, 0.86, 0.91, 0.94, 0.97), negative stel- 
lar masses were obtained if a (dyn) (i?, v?) = 0. The condition 
M(d, cp) > gives an upper limit for a' dyn ^ (R, cp) < at 
these phases and shows clearly the limitations of treating 
an RR atmosphere in QSAA. In these phases the character- 
istic value \q\ w 1. The sharp peaks q(R, cp = 0.5) = 38.8 and 
q(R, cp = 0.65) = —13.8 are produced by dv/dt « —vdv/dr, 

1. e. the atmosphere is almost completely free of acceleration. 
q(R,cp = 0.05) = 11.4 is a result from the poor representa- 
tion of the density stratification by (0), i.e. the dynamical 



term is large in g c — g s 



,(dyn) 



= (44.7 - 8.7 - 33.1) re- 



producing a (dyn) (R, cp) ss ll(dv/dt + vdv/dr). This is the 
start of the descending branch in the light curve, when the 
atmosphere starts expanding rapidly. 

The averaged surface gravity for the whole pulsation 
cycle is (log<7 e ) = 2.73, which is in excellent agreement with 
(log^ge) = 2. 69, 2.72 from uvby and UBV(RI)c photometry 



Sicgcl 19821. IlTu fc J anes 1990j), respectively. The averaged 
effective temperature of (T e ) = 67 78 K is signif i cantly higher 
than (T e ) = 6433 6400, 6490 K (|Siegellll982l . iLiu fc Janes! 
Il990l , |Barczal l2003). The minimal, averaged, and maximal 
angular radii are (1.519,1.669,1.801) x 10" 10 rad corre- 
sponding to .Rmin = (4.46, 4.90, 5.29)i?Q , respectively. The 
following equilibrium luminosi ty and effective temperature 
l|Carnev. Strom fc Jones |[l992l ') give the position of SU Dra 
in a theoretical HRD: L eq = 4-wod 2 {if (cp)T*(cp)) = (45.9 ± 
9.3)L Q , T cq = {L c j4namd] 2 }^ 4 = {i9 2 T , c 4 ) 1//4 (#)~ 1 / 2 = 
(6813 ± 20)K. The magnitude averaged absolute brightness 
is (My) = +0.68 ± .23 mag. 



4 DISCUSSION 

From the point of view of radiative transfer, th e plane- 
paral lel approximation of the model atmospheres dKurucj 
Il997l ) could be applied beyond doubt, because h^ 1 < 0.037? 
holds for the whole pulsation cycle of an RR star. In favour of 
applying QSAA, semi-quantitative arguments were that the 
temperature changes are very slow even in a high amplitude 
RR star like SU Dra: 5K/2500s < \dT c /dt\ < 1500K/5000s, 
while the chara cteristic timescale of the radiative processes 
is be low 200s (|Oke. Giver fc Searlel 1 19621 . iBuonaura et al.l 
1985). The effect of the variable effective gravity can be 
characterized by our Condition II, which provides infor- 
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mation whether QSAA may be assumed. The limits are 
4.5cms _2 /2500s < \dg e /dt\ < 14200cms _2 /2500s. The up- 
per limit comes from the rising branch of the light curve, 
when the main shock hits the atmosphere. 

In general, the photometric input of the present method 
is identical with that of the BW method. In order to obtain 
the fundamental parameters, radial velocity data and their 
problematic conversion to pulsation velocities are not nec- 
essary. On the other hand, $ and ho must be differentiated 
numerically, differential quotients are sensitive to the non- 
validity of QSAA. 

Our photometric and hydrodynamic considerations re- 
vealed empirical quantitative conditions to find the phase 
int ervals o f the p ulsation when static model atmospheres 
of H urucj l| 19971 ) are satisfactory to derive the variable 
and non-variable physical parameters of the pulsating at- 
mosphere. Outside these intervals, dynamical model atmo- 
spheres are necessary to refine the parameters from QSAA, 
which is beyond the scope of the present paper. The funda- 
mental parameters d, M, [M], E(B — V) were determined 
using photometric quantities only in phases when the QSAA 
is a good approximation (i.e. both Conditions I and II were 
satisfied). The values obtained from averaging over the entire 
pulsation cycle can be considered as a first approximation 
only, because QSAA was assumed in all phases regardless of 
it being a good or poor approximation. The large error of 
L eq and (My) originates from Ad/d w .1 of our best value 
[*] in Table [2] 

To give an impression on the accuracy of inverting the 
UBV(RI)c photometry to physical parameters, the com- 
parison star BD +67 708 was used, becaus e its colours are 
similar to those of SU Dra (|Barczal [20o3 . Table 3). The 
results are [M] = -0.77 ± .03, E(B - V) = .000, = 
(1.913±.002)xl0~ 10 rad, \ogg = 3.59±.01, T c = 7505±5 K. 
The errors are roughly in the same order of magnitude as 
those of SU Dra when Conditions I and II are satisfied. 

4.1 Kinematics of the atmosphere 

Figs. [3J;, d demonstrate that significant corrections must 
be added to 7?, R if the true pulsation velocity and ac- 
celeration are required at < r < 1. Moreover, the tri- 
angles show another correction — 4.6kms _1 < v(R,ip) = 
h^dv/dr < S.Skms -1 if 0.1 < p < 0.6, i.e. v(r,tp) = 
v (R, tp) — v(R, tp)(R — r)ho(R, tp) H — • is the velocity profile. 
In other words, even in the shock free phases, considerable 
phase-dependent velocity and acceleration gradients exist in 
the layers R — hg 1 < r < R of the line formation. The prob- 
lem of v ^ is present in all phases, i.e. even in the shock 
free intervals, whi ch are used in modern BW analyses (e.g. 
iLiu fc JanesHl990T ). 

The non-uniform motion of the outermost layers intro- 
duces uncertainty when pulsation velocities are determined, 
because the centre of mass velocity n 7 must be subtracted 
from the observed radial velocities. A rece nt exposition of 
the p roblem for Cepheid stars is given in iNardetto et al.l 
l|2009i ). By definition refer to < t « 1, while the 

variable component of the radial velocity is an average of 
velocities (e.g. (|9])) over the layers R — < r < R i.e. < 
t < 0.3. The effect on « 7 and d has not yet been studied at 
all. Nevertheless, the importance is obvious, sinc e an error 
of 1 kms -1 in u 7 results in an error Ad/d ~ 0.1 (|Gautschvl 



Il987h . There is a considerable uncertainty of v 1 in the litera- 
ture, suggesting an error of d as large as a factor of 2. It suf- 
fices to ment ion that for SU Dra t > 7 = — 161, — 166.9kms~ 
are gi ven by lOke. Giver fc Searld il962h and ILiu fc Janesl 
( 1990) from high dispersion spectra and spectral masking 
method CORAVEL, respectively. 

A wavy fine structure in the variation of R is clearly 
seen in Figs. [T] |3p. Standstills are at tp w 0.25, 0.55, 0.78, 
when the outward motion is reversed, while the outward mo- 
tion starts after short standstills at (p » 0.45, 0.74, 0.94-0.98. 
They are without doubt real, because the beginnings of 
the outward motion are connected with significant max- 
ima in T e (<p), log g e (tp)- The well k nown bump and hump 
at ip « 0.74,0.94-0.98 (|Smithlll995l ) in the light curve are 
caused by the precursor and main shocks, respectively. At 
ip ~ 0.45, a small change is observ able in the slope of the 
light curve l|Barcza fc Benko]|2009r i. Here, the existence of 
a shock is a new finding. It is a pre-precursor shock; we 
propose the designation jump for it. By integrating radial 
velocities, authors tend t o smooth out the m entioned fine 
structure of motions (e.g. ILiu fc J anes 1990 in the case of 
SW And). This practice seems to be unjustified. 

In the interval <p « 0.82-0.90, the atmosphere is roughly 
in standstill, the brightness starts rising, a small depression 
of &(<p), \ogg e (tp) is visible at ip w 0.84-0.86, but T a (tp) is 
monotonic. It is not clear whether the small undulation of 
$(<p) around (1.571-1.579) x 10~ 10 is real or not, because the 
maximal A log g c — 0.116 was found just at tp = 0.88. 

4.2 T c scale, mass, distance, absolute brightness 

Surprisingly, in comparison with previous studies, consider- 
able differences were found only in (T c ), T oq . This is due to 
the fact that the interpolated T c (tp) depends on log g e (<p) 
and that the information from a five colour photometry was 
used in a more complex manner: averaged value from 30 
colour index pairs was utilized. The photometry covers the 
whole spectrum between 350 and 1000 nm. The use of one 
colour in dex only, e.g. V — K solely "because of its apparent 
merits" (ILiu fc Janes!ll990h . can result in systematic error 
of T e (tp). It can explain the « 350 K difference, because in 
that previous work QSAA was assumed in all phases regard- 
less of violating both Conditions I and II. An inspection of 
the functions {T C W (log q e , CIi , Clg , \M ] , E(B - V ) ) } i=1 , 2 de- 
rived from the tables |Kuruczlll99"7h shows that, if merely 
one colour index is to be used, the optimal choice for deter- 
mining T e would be Rc — Ic- This is because Rc — Ic is 
almost independent of log g c for the actual values of [M] and 
E(B — V) of SU Dra. However, in phases violating Condition 
I, the sole use of Rc — Ic would also introduce a systematic 
error, similarly to the use of V — K . 

There is a remarkable decrease of Ad = 200 — > 67pc, 
AM = 0.10 — > O.OSMq if our compressible QSAA is substi- 
tuted for UAA, i.e. more physical input is used in the frame 
of a one-dimensional model in space. Since (|10|l was derived 
from the dynamic equation ([3}, the mass determination is 
more accurate from it while the distance is more uncertain. 
This is reflected in the shape of the curves in Fig. [3^,. The 
large formal error of d originates from the non-separability 
of errors in #, ho, dho/dt, d 2 h /dt 2 , and a (dyn) (r,t). 

A first attempt to derive d and M of an RR st ar from 
the over-simplified version |T} of ((3| was described bv lBarczal 
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The present results for the distance and mass of SU 
Dra differ only slightly from the previous d = (647 ± 16)pc, 
M — (0.66 ± .O3)A^0 . The most probable reason of the very 
good coincidence is the existence of the phase island with 
the good QSAA at ip « 0.93, just when the atmosphere is in 
standstill, consequently JTJ is a good approximation because 
of v x 0. 

The values d = 640 pc, M = 0.47X Q given by 
iLiu fc Janes! (|l990T ) are very close to the present ones. How- 
ever, some caution is appropriate because of the underesti- 
mated uncertainties in their derivation. Their M originates 
from using fTJ, i.e. C — 0, UAA corresponding to line 3 in 
Table [5] Furthermore, because of the uncertain value of v 1 , 
Ad/d < 0.6 is well possible. 

In the BW method, the propagation of the error Ai> 7 = 
1 kms -1 can be estimated from ([2]) as follows. Typical ra- 
dius changes of an RR star are 7? ma x — -Rmm ~ 5 x 10 km 
within P/2 ~ half day. The error of the radius change is 
AR « Av 1 P/2 « 5 x 10 4 km, i.e. AR/(R max - R mitx ) « 0.1. 
By d — (i?max - Rmia ± A7?)/(i? max - i? m in), the final er- 
ror will be Ad/d « 0.1. This considerably ex ceeds the er- 
ror o riginating from the projection factor (e.g. ILiu fc Janes! 
1990) converting the observed radial velocity to pulsation 
velocity. Furthermore, by ((2| and UJ, the error Ad/d w 0.1 
propagates to an error AM/M w 0.2. 

The lower limit of the magnitude averaged visual ab- 
solute brightness is My > +0.48 mag if d < 724 pc, while 
Mv = +0.26, +0.01 mag belong to the improbable values 
d = 800 and 900 pc, respectively. 



5 CONCLUSIONS 

Observed colours and magnitudes of a spherically pulsating 
star have been compared with those of static Kurucz model 
atmospheres to determine fundamental parameters of the 
star in the frame of the quasi-static atmosphere approxima- 
tion. Photometric and hydrodynamic conditions have been 
formulated for the validity of the quasi-static atmosphere 
approximation in spherically pulsating stars. 

(1) The quasi-static atmosphere approximation has been 
generalized for a non-uniform, compressible atmosphere 
with radial velocity gradient. This is a step forward be- 
cause the hitherto available uniform atmosphere approxi- 
mation described the motion of the atmosphere by effective 
gravity and differentiating the sole parameter radius R of 
optical depth zero with respect to time. 

(2) Our combined photometric and hydrodynamic method 
uses the variation of effective gravity, angular radius, veloc- 
ity, acceleration in the Euler equation. The complete input 
of the Euler equation has been derived from photometry 
and theoretical model atmospheres. Spectroscopic and ra- 
dial velocity observations were used neither explicitly nor 
implicitly. This is a definite advantage in comparison with 
the Baade-Wesselink method, because less observational ef- 
forts are needed to get the fundamental parameters and the 
problematic conversion of observed radial velocities to pul- 
sation velocities is not necessary. 

(3) Concerning the interpretation of multicolour photom- 
etry, the inputs are identical with those of the Baade- 
Wesselink method. The present refinements are the quan- 



titative conditions whether photometry can or cannot be 
interpreted by static model atmospheres. 

(4) Firstly, the phases have been selected by the photomet- 
ric conditions when the quasi-static atmosphere approxima- 
tion is valid. Secondly, in a number of these phases, the laws 
of mass and momentum conservation have been applied in 
Euler formalism of hydrodynamics to determine mass and 
distance of the star from the motion of the atmospheric lay- 
ers in the neighbourhood of zero optical depth. Afterwards, 
it was checked whether the hydrodynamic condition of the 
quasi-static atmosphere approximation was satisfied in the 
used phases, i.e. phases that satisfied the photometric con- 
dition but violated the hydrodynamic condition had to be 
excluded. Atmospheric dynamical mass seems to be an ap- 
propriate term to indicate that the mass has been derived by 
a method which is completely diffe rent from a m ass derived 
by pulsation or evolution theories (Smith 1995). 

(5) As a by-product, a variation procedure has been given 
for estimating atmospheric metallicity and interstellar red- 
dening toward a star from UBV(RI)c photometry. This 
method was successfully applied for the non-variable com- 
parison star BD +67 708, giving [M] = -0.77 ± .03, E(B - 
V) = .00, = (1.913 + .002) x 10" 10 rad, log g = 3.59 ±.01, 
T c = 7505 ± 5 K. 

(6) Using the UBV(RI)c photometry of the high ampli- 
tude RRab star SU Dra, the following fundamental param- 
eters have been found: 

[M] = -1.60 ± .10, 

E{B-V) = 0.015 ±.010, 

d = (663 ± 67)pc, 

Rutin — 4.46J?0, -Rmax = 5.29-Rq , 

M = (0.68 ± M)M Q , 
L cq = (45.9±9.3)L Q , 
T cq = (6813 ± 20)K. 

L oq and T eq are approximate values, since they originate 
from averaging over the whole pulsation cycle containing 
phases in which the quasi-static atmosphere approximation 
is merely a first approximation. 

(7) The internal motions of the atmosphere with respect 
to zero optical depth have been found to be significant in 
comparison with velocities and accelerations derived from 
the uniform atmosphere approximation. From the inter- 
nal motions of the atmosphere, some constraints have been 
sketched for converting observed radial velocities to pulsa- 
tion and centre of mass velocities. Estimations have been 
given for the error propagation in a Baade-Wesselink anal- 
ysis. 
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